Free-energy landscape of nucleation with an intermediate metastable phase studied 

using capillarity approximation 
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Capillarity approximation is used to study the free-energy landscape of nucleation when an in- 
termediate metastable phase exists. The critical nucleus that corresponds to the saddle point of 
the free-energy landscape as well as the whole free-energy landscape can be studied using this cap- 
illarity approximation, and various scenarios of nucleation and growth can be elucidated. In this 
study we consider a model in which a stable solid phase nucleates within a metastable vapor phase 
when an intermediate metastable liquid phase exists. We predict that a composite critical nucleus 
that consists of a solid core and a liquid wetting layer as well as pure liquid and pure solid critical 
nuclei can exist depending not only on the supersaturation of the liquid phase relative to that of 
the vapor phase but also on the wetting behavior of the liquid surrounding the solid. The existence 
of liquid critical nucleus indicates that the phase transformation from metastable vapor to stable 
solid occurs via the intermediate metastable liquid phase, which is quite similar to the scenario of 
nucleation observed in proteins and colloidal systems. By studying the minimum-free-energy path 
on the free-energy landscape, we can study the evolution of the composition of solid and liquid 
within nuclei not limited to the critical nucleus. 

PACS numbers: 64.60.-i, 64.60.Q-, 82.60.Nh 
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I. INTRODUCTION 

A phase transformation involves nucleation and 
growth. When an intermediate metastable phase exists 
between the initial metastable phase and the final sta- 
ble phase, the scenario of phase transformation becomes 
complex and has attracted much interest from more than 
a decade agoi. Recently, renewed interest has emerged 
not only in the field of traditional metallurgy^ but also 
in the field of soft-condensed-matter physics of proteins 
and colloids?"--. Even for the simple Lennard- Jones sys- 
tem below the triple point, it has recently been suggested 
— that the intermediate metastable liquid plays a crucial 
role in the vapor to solid phase transformation. 

The direct microscopic computer simulation of a phase 
transformation using molecular dynamics or the Monte 
Carlo method is possible^ but is still a difficult task. To 
avoid the requirement of huge computational resources 
and to obtain a qualitative (coarse-grained) picture of 
the kinetics of a phase transformation, a mesoscopic ap- 
proach called the phase-field model, which is based on 
density functional theory^—, has been frequently used. 
Using this phase-field model for the nonconserved or- 
der parameter, Bcchhoefer et al»^i and Celestini and ten 
Boschi^ discovered the formation of a finite layer of an 
intermediate metastable phase at the growing front of 
the stable phase. Their results were also confirmed for 
a different class of free-energy landscape^. Later, these 
works were further extended to the conserved order pa- 
rameter-^. These previous works, however, focused on 
growtbii"— rather than nucleation. Usually the exis- 
tence of a critical nucleus has been assumed from the 
outset of the simulation^. 

It has been customary to assume two-step nucleation 
when there is an intermediate metastable phased—. 



Imagine that nucleation proceeds via two successive pro- 
cesses: the nucleation of a metastable liquid nucleus 
within a vapor phase, and the subsequent nucleation of 
a solid phase within the metastable liquid nucleus. Fig- 
ure [3a ) shows a schematic diagram of the free energy G 
of nucleation along a fictitious one-dimensional reaction 
coordinate when an intermediate metastable phase ex- 
ists where a metastable vapor phase (V) transforms into 
a stable solid phase (S) via an intermediate metastable 
liquid phase (L). This diagram implicitly assumes that 
the reaction coordinate is one-dimensional. Also, the nu- 
cleus of the stable S phase is assumed to grow within the 
nucleus of the metastable L phase (Fig. D^b)), implying 
successive nucleation and growth. 




FIG. 1: (a) Schematic of the free energy G of a three-phase 
system along a fictitious reaction coordinate, (b) Composite 
nucleus consisting of a stable solid core surrounded by an 
intermediate metastable liquid layer in a metastable vapor 
phase. 



Kashchiev and coworker a 15 ' 16 have assumed the above 
scenario and considered both nucleation and growth on 
the same footing within the framework of Kolmogorov- 
Johnson-Mehl-Avrami (KJMA) kinetics^—. They as- 
sumed that these two nucleation processes occur suc- 
cessively (Fig. QJb)) and derived a formula for the time 
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dependence of the volume fraction of the stable phase. 
They found, in particular, that the nucleation of the sta- 
ble phase can be delayed by the existence of an interme- 
diate metastable phase. Their result also suggests the ap- 
pearance of the long-lived metastable intermediate phase 
predicted by the phase-held model^ ' 11 ' 12 ' 14 . Similarly, 
Valencia and Lipowsk y 17 ' 18 derived a formula for the nu- 
cleation rate assuming the intermediate metastable phase 
using the double-barrier model by extending the theory 
of Kramers on the stochastic process of nucleatio n 23 ' 24 . 
Their formula also suggests that the nucleation rate of 
the stable phase is lower if the intermediate metastable 
phase exists. On the other hand, Nicolis and Nicolis 1 ^ 
showed that the nucleation rate of the stable phase can 
be enhanced by the existence of the metastable phase by 
using the theory of Kramer o 23 ' 24 . A similar enhancement 
of the nucleation rate by the existence of an intermediate 
metastable phase was directly observed by a numerical 
experiment by Sear— based on the lattice model . In 
these theoretical studies 1 ^— , however, it was implicitly 
assumed that the nucleation reaction coordinate is one- 
dimensional (Fig. da)) and that there are two distinct 
nucleation barriers (a double barrier). 

In fact, nucleation proceeds on the free-energy surface 
in multidimensional phase space, and the nucleation bar- 
rier corresponds to the saddle point on the multidimen- 
sional surfac o 26 ' 27 . It is, therefore, not apparent if two 
distinct saddle points really exist in the free-energy land- 
scape as assumed by previous authors 1 ^— . Incidentally, 
our problem resembles that of the vapor phase nucleation 
of a binary droplet and, therefore, resembles that of the 
deliquescence 2 ^—. In our problem, however, two com- 
ponents, the liquid and the solid, come from the same 
metastable vapor phase. On the other hand, the solute 
and the solvent in the deliquescence come separately from 
the pre-existing solid core and vapor phase respectively. 

In this study, we use the classical capillarity approx- 
imation based on the assumption of uniform density to 
reduce the multidimensional free-energy surface to a sur- 
face in a finite parameter space and study the nucleation 
pathway when an intermediate metastable phase is in- 
volved. We use a crude capillarity approximation that 
is also the basis of classical nucleation theory (CNT) as 
it allows us to study the critical nucleus as well as the 
whole process of nucleation of a single nucleus. 



II. CAPILLARITY THEORY OF NUCLEATION 
WITH AN INTERMEDIATE METASTABLE 
PHASE 

We hrst consider a composite nucleus consisting of a 
solid core with radius r embedded in a liquid nucleus 
of radius R (Fig. [2(b)). This problem is akin to sur- 
face melting or wetting 3 - 1 - - — , where the thickness of the 
wetting layer is given by R — r. In fact, our problem is 
inverse melting, where a stable substrate appears after 
the appearance of the wetting layer. 



Within the capillarity approximation, the free energy 
required to form a liquid nucleus with radius R from a 
metastable vapor phase is given by 

AG LV = A r (ML-Mv)+47ri? 2 a LV , (1) 

where R is the radius of the liquid nucleus, N is the num- 
ber of liquid molecules within radius R, ctlv is the liquid- 
vapor surface tension, and /il and /iy are the chemical 
potentials of the liquid and vapor phases, respectively. 
Similarly, the free energy required to form a solid nu- 
cleus of radius r within a metastable liquid nucleus of 
radius R (Fig.QTb)) is given by 

AG SL = (iV-7V 1 )(/x s -/x L ) + 47rrVsL ! (2) 

where N — Ni is the number of solid molecules, Ni is 
the number of liquid molecules in the liquid layer with 
width R — r, /is is the chemical potential of the stable 
solid phase, and <tsl is the solid-liquid surface tension. 

The total free energy AC? required to form a composite 
nucleus with a solid core of radius r surrounded by a 
liquid layer of width R — r is given by 

AG = AGlv + AG SL + AGslv, (3) 

where 

AGslv = 4 7 ri? 2 S , e-( fl - r )/« (4) 

is the correction term due to the short-range inter- 
action between the liquid-vapor and solid-liquid inter- 
face s 31 ' 33 ' 34 which corresponds to the disjoining pressure 
in the theory of deliqescenc o 29 ' 30 , and 

S = osv - olv - c S l (5) 

is the spreading parameter—, where asv is the solid- 
vapor surface tension, and £ specifies the range of in- 
teractio n 31 i 33 ' 34 . We have used the crudest possible ap- 
proximation to the short-range interactions for wetting 
in Eq. since we are most interested in the global pic- 
ture of nucleation, although a more rchned theory for 
short-range wetting is available i 3 - 7 - 
From Young's equation 

csv = osl + ctlv cos 6, (6) 

where 9 is the contact angl o 35 ' 36 , we have 

S = cr L v (cos 6»-l). (7) 

Therefore, the complete wetting of the solid by the liquid 
is realized when S > 0, otherwise the solid is incompletely 
wet by the liquid layer. Even though Eq. ([7]) predicts that 
the complete wetting with 9 = implies S = 0, it is well 
recognized that the surface free energy asv, fSL and ctlv 
will change 3 -^ even in the complete wetting regime with 
9 = such that S becomes positive (S > 0). Physically, 
Eq. represents the interaction between the solid- liquid 
and liquid- vapor interfaces, which is repulsive (S > 0) 
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when the complete wetting condition is satisfied as these 
two interfaces repel each other so that the liquid wetting 
layer intervenes between the solid and vapor phases. The 
condition S > is also known as the condition of surface 
melting 3 ^. 

From Eqs. {TJ to (gj), AG in Eq. © becomes 

AG = N (ml - Mv) + (N - JVi) (ms - Ml) 

+ 47ri? 2 (a LV + (r 2 /i? 2 ) + fl-e-^--)^) ,(8) 

which leads to the free energy of a solid nucleus directly 
nucleated from the vapor when r = R and N\ = 0: 



AG = AGsv =N(ns- Mv) + 47ri? 2 o- S v- 



(9) 



By using N = AttR 3 /3v m and maximizing this free energy 
by solving 9 AGsv /dR = 0, we obtain the free-energy 
barrier AGg V and the critical radius i?g V of the solid 
critical nucleus: 



AG 



sv 



167T <g|y 

3 A/i§ v ' 



R 



sv 



A^sv 



(10) 



where v m is the molecular volume and A/isv = Mv~Ms > 
0. Similarly, the free energy of a liquid nucleus is obtained 
when r = and JVi = iV: 

AG = = AGlv = N (ml - Mv) + 4?ri? 2 (a LV + Se-*^ 

~ iV(^ L - J u v ) + 47ri? 2 ( j L v, (11) 

where the approximation of the second line is valid pro- 
vided £ <C R- The free-energy barrier AG£ V an d criti- 
cal radius i? LV are given by formulae similar to those in 
Eq. (TIT)]) obtained by changing the suffix from "SV" to 
" LV" and using the chemical potential difference Amlv = 
Mv — Ml • The molecular volume v m is assumed to be the 
same in the solid and liquid nuclei. 

However, instead of using the above expressions for 
the free energy, we introduce the following scaled free 
energies: 

g = AG/AG* SV , g LV = AG LV /AG* SV , 
g SL = AG SL /AG SV , 5 sv = AGsv/AG^ v , (12) 
.9slv = AGsLv/AGg V , 

and express these energies in terms of two parameters: 
the scaled radius x of the liquid nucleus, defined by 



= R/R. 



SV; 



(13) 



and the proportion t of the solid radius r relative to the 
liquid radius i?, 



t = r/R. 



(14) 



By varying t in the range < t < 1, we can study a 
composite nucleus composed of a solid core surrounded 
by a liquid wetting layer. The nucleus is all solid when 
t = 1 and all liquid when t = 0. 



The above free energies as functions of the two param- 
eters (x, t) are given by 



g(x,t) 



-2Sx 3 - 
Zx 2 ( 8- 



2(1 



at 2 



5) (xtf, (15) 
(1 - {a + B^e-^ 1 -^/^ 



and 

ffsv 

.9LV 

5s L 
5s LV 



-2a; 3 + 3a; 2 



g (x, 0) = -2Sx 3 + 3a; 2 (/? + (1 - (a + B)) e~ x 'A , 
-2 (1 - 5) (xtf + 3a (xt) 2 , (16) 
3a; 2 (1- (a + 8)) L-*<l-*)h 



,-x/t 



where we have included the exponential correction of 
AGslv (Eq. I@J)) in the definition of <?lv (the last term 
of Eq. pi[) ). and introduced the material parameters 



a = osl/osv, /3 = o- L v/csv, 
S = A^ L v/A/isv, r = £/-Rg V - 



(17) 



The total free energy of the composite nucleus 
[Eq. (|15p] is rewritten as 



9 = 9lv + gsh + 5SLV- 
The spreading parameter S is given by 
S/vsv = l-(a+B), 



(18) 



(19) 



using the reduced surface tensions a and 8. From the 
Dupre equation^, a SL « cr S v + olv - V^svOlv = 

(v^sv — x/f lv) ' we have a « (l — \/p) ■ Since the 
density of the solid phase is close to that of the liq- 
uid phase except near the critical point, we expect that 
ctsv ~ °lv, ft ~ 1) and a <C 1. Incomplete wetting 
(S < 0) is realized when a + 8 > 1. 

The critical radius Xg V and activation energy gjy of 
a solid nucleus that is directly nucleated from vapor are 
given by Eq. (|10[) and are written as 



1. 



9sv 



1. 



(20) 



Similarly, formulae for a metastable liquid nucleus in va- 
por can be obtained from dgi_y/dx = and are approxi- 
mately given by 



B 3 



5' 5lv "^- 



(21) 



On the other hand, those for a solid nucleus that is nu- 
cleated from a metastable liquid phase can be obtained 
from dgsh/dt = and are given by 



t 



SL 



a;(l-<5) ; 



9sl 



(1-6) 



2 ■ 



(22) 



If B 3 /5 2 > 1, the initial critical nucleus will be mostly 
solid as <7£ v > <?gy This scenario is expected when 8 ~ 
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1 > 8 and the metastable liquid phase is closer to the 
metastable vapor phase than the stable solid phase (8 = 
Afius / Afisv < !)• Since the liquid phase is less stable, 
the solid phase will be directly nucleated from the vapor 
phase. 

On the other hand, when f3 3 /S 2 < 1 the initial critical 
nucleus will be liquid. In this case, the metastable liquid 
phase is closer to the stable solid phase (8 w 1). A stable 
solid nucleus will grow within the metastable liquid ma- 
trix, and two-step nucleationi^— with double barriers is 
expected because the most stable phase is the solid. 

Therefore, by increasing the relative supersaturation 8 
of the liquid phase, one may expect the bifurcation from 
a solidlike critical nucleus surrounded by a thin wetting 
layer of liquid for small S to a liquidlike critical nucleus for 
large 8. This solidlike to liquidlikc change of the charac- 
ter of the critical nucleus was theoretically predicted by 
Granasy and Oxtoby-2 using the triple-parabola model, 
and is implied by Fig. 4 of Talanqucr and Oxtoby^ using 
density functional theory. However, they only considered 
the critical nucleus at the saddle point of the free-energy 
landscape because density functional theory can only be 
used to study the critical nucleus that corresponds to the 
stationary state at the saddle point on the free-energy 
landscape. The evolution of the nucleus and its compo- 
sition other than that of the critical nucleus can only be 
determined through the examination of the free-energy 
landscape. To study the nucleation scenario in the free- 
energy landscape qualitatively, we apply our capillarity 
theory in the next section. 




(b)<5 = 0.5 



III. RESULTS AND DISCUSSION 
A. Incomplete wetting (a + /3 > 1) 

In this case, the metastable liquid may not wet the 
solid nucleus. Therefore, it will be unfavorable for the 
solid phase to nucleate within the liquid phase. Thus, 
the composite nucleus is not expected to appear. Fig- 
ures [DJa) and (b) show contour plots of the free-energy 
landscape of g(x,t) ((a) 8 = —0.2, (b) 8 = 0.50) in the 
x — t plane when the chemical potential /il of the liquid is 
closer to that (/Liy) of the vapor than the solid (/lis)- The 
nucleation pathway starts at any point along x — and 
ends at x = oo and t — as the solid phase is the sta- 
ble phase. We used the reduced surface tensions a = 0.3 
and /3 = 0.8, which correspond to the incomplete wetting 
condition S > 0. The parameter t is fixed to r = 0.2. 

Note that the pathway along the t = axis is the liquid 
axis, which corresponds to the nucleation of a metastable 
liquid nucleus, and that along the t = 1 axis is the solid 
axis, which corresponds to the direct nucleation of a solid 
nucleus from the vapor. A negative 8 means that the liq- 
uid phase is less stable than the vapor phase. In both 
Fig. [5{a) and Fig. H^b) we can locate the saddle point 
at x ~ 1 and t = 1, indicated by red points on the 
minimum-free-energy paths (MFEPs) indicated by the 



FIG. 2: Contour plots of the free-energy landscape of g(x,t) 
when the relative supersaturation is low: (a) S — —0.2 and 
(b) 8 = 0.5. The metastable liquid phase cannot appear and 
a solid nucleus directly grows within the metastable vapor 
phase. The other parameters used are a = 0.3, /3 = 0.8, 
and r = 0.2. The red solid line along t = (the solid axis) 
indicates the minimum-free-energy path (MFEP) of the nu- 
cleation process. The saddle points are indicated by the red 
point on each path. 



red solid lines. This means that the solid phase nucle- 
ates directly from the vapor without passing through the 
liquid phase. Since the metastable liquid phase is closer 
to the vapor phase and is less stable, the metastable liq- 
uid phase cannot appear. Instead, a solid nucleus directly 
nucleates from the vapor phase as expected from the dis- 
cussion in the last section. This solid nucleus does not 
accompany the liquid wetting layer as the liquid does not 
wet the solid. 

Figures [Ha) and (b) show contour plots of the free- 
energy landscape when the relative supersaturation 8 of 
the liquid phase is higher (8 = 0.72 and 8 = 0.80). Pa- 
rameters a, j3, and r are the same as those in Fig. [5] In 
Fig.(3Ja), we can clearly observe the existence of two sad- 
dle points on the solid axis (t = 1) at (x,t) = (1.0,1.0) 
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and the liquid axis (t = 0) at (x,t) = (1.1,0.0) indicated 
by red points on the MFEPs. The free-energy barrier at 
these two saddle points is the same (g* = 1.0). Therefore, 
two nucleation pathways (MFEPs), indicated by the two 
solid lines on the liquid (t = 0) and solid (t = 1) axes, 
coexist. 




0.5 1 1.5 

x 

(a)<5 = 0.72 




FIG. 3: Contour plots of the free-energy landscape of g(x,t) 
for the nucleation of the composite nucleus for the moderate 
and high relative supersaturations: (a) S — 0.72 and (b) <5 = 
0.80 when a = 0.3, /3 = 0.8, and r = 0.2. The red solid 
curve along the solid axis (t — 0) indicates the MFEP of 
the nucleation process. The red broken line along the liquid 
axis (t = 1) indicates another MFEP. The saddle points are 
indicated by the two points on the paths. 



Of course, the MFEP along the liquid axis will eventu- 
ally merge into the MFEP along the solid axis far from 
the saddle point as the free energy of the solid phase at 
(x — > oo,i = 1) should be lower than that of the liquid 
phase at (x — > oo,i = 0). 

As the relative supersaturation S is further increased 
and the metastable liquid phase becomes more stable, 
the saddle point at (x,t) = (1.0,0.0) on the liquid axis 
becomes lower (g* = 0.80) than that (g* = 1.0) at 
(x,t) = (1.0,1.0) on the solid axis (Fig. E^b)). Fig- 



ure 2Ja) shows cross sections of the free-energy surface 
along the solid axis (t = 1) and liquid axis (t = 0) when 
S = 0.72. When S becomes larger and the liquid phase 
becomes as stable as the solid phase, the liquid nucleus 
along the red solid line on the liquid axis becomes more 
probable than that along the red broken line on the solid 
axis (Fig. [3Jb) ) . Again, the MFEP along the liquid axis 
will eventually turn toward that along the stable solid 
phase at (x — > oo, t = 0) far from the saddle point 
as shown in Fig. 0|b). Therefore, the nucleation path- 
way from metastable vapor to stable solid through the 
intermediate metastable liquid state is free-energetically 
easier than that from vapor to solid directly. Similar re- 
sults that the nucleation pathway with intermediate state 
has a lower free energy are obtained by ten Wolde and 
Frcnkei£ using Monte Carlo simulation and by Lutsko 
and Nicolis^i using density functional theory. 




(a)Cross sections along t = and t = 1 
(<5 = 0.80) 




(b)MFEP of supercritical nucleus (<5 = 0.80) 



FIG. 4: (a) Cross sections of the free-energy surface along 
the solid axis (t = 1) and liquid axis (t = 0) when the relative 
supersaturation is high (S = 0.80). Clearly the saddle point 
(maximum) of the free-energy curve along the liquid axis is 
lower than that along the solid axis, (b) Contour plot of 
the free-energy landscape of g(x,t) and the MFEP (red solid 
curve) for the supercritical nucleus. The MFEP eventually 
converges to the upper right corner of the pure solid phase. 



There are now two pathways of nucleation. One is the 
direct nucleation of the solid phase from the metastable 
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vapor. The other is the indirect nuclcation of the 
metastable liquid phase from the vapor, which will even- 
tually transform into the solid phase without crossing the 
barrier as shown in Fig. Hth). Therefore, there are two 
nucleation rates. One is the vapor to liquid nucleation 
rate, 



Jlv = A h y exp 



kT 



(23) 



which expresses the nucleation rate of a liquid nucleus, 
where AG LV is the free-energy barrier for the liquid nu- 
cleus, given by Eq. (|2ip . and kT is the temperature. Alv 
is the so-called prccxponcntial factor which is the kinetic 
factor of molecules attaching to and detaching from a 
nucleus. The other is the vapor to solid nucleation rate, 
given by 



Jsv = A sv exp 



AG 



sv 



kT 



(24) 



where AGg V and Asv are the free-energy barrier and the 
preexponential factor similar to those in Eq. (|2~3l) . 

Previously, it has been assumed that nuclcation oc- 
curs on a one-dimensional energy landscape through the 
successive nucleation of the liquid from vapor with nucle- 
ation rate Jlv followed by liquid to solid nucleation with 
rate J S l, given by 



Jsl = ^sl exp 



AG 



SL 



kT 



(25) 



where AGg L and ^4sl are the free-energy barrier and the 
preexponential factor similar to those in Eq. (|23|) . Using 
the kinetic theory of nucleation, where the total balance 
of attachment and detachment of molecules is consid- 
ered^, a formula similar to the conductance of two reg- 
isters connected in series applies approximately for the 
total nuclcation rate J—: 

J Jlv Jsl ^ ^ 

where only the contribution around the two saddle points 
for JgL and Jlv are retaine d 17 ' 18 . Therefore, the total 
nucleation rate J is slower (lower) than Jlv- However, 
in our model two nucleation processes occur in paral- 
lel. Also, liquid to solid nucleation occurs after vapor 
to liquid nucleation without crossing the energy barrier. 
Therefore, a formula similar to the conductance of two 
registors connected in parallel follows: 



J PS J 



LV 



Jsv j 



(27) 



where the interference of two channels shown in Figs. [3th) 
and EJb) as the solid and the broken line is neglected. 
Then, the total nuclcation rate J is faster (higher) than 
Jlv- The nucleation can be enhanced because not only 
a direct vapor to solid nucleation channel exists but also 
the vapor to liquid nucleation is followed by a barrierless 
liquid to solid transformation. A similar but slightly dif- 
ferent explanation of the enhancement of the nuclcation 
rate due to the presence of a metastable phase near the 
spinodal has been proposed^. 



B. Complete wetting (a + f3 < 1) 

The critical nucleus is either solid or liquid for the in- 
complete wetting case because a composite nucleus that 
consists of a solid nucleus surrounded by a metastable 
liquid wetting layer is energetically unfavorable. In con- 
trast, such a composite nucleus is expected to occur for 
the complete wetting case. 




0.5 1 1.5 



(a)<5 = 0.3 




(b)<5 = 0.5 



FIG. 5: Contour plots of the free-energy landscape of g{x,t) 
when the relative supersaturation is low: (a) S = 0.3 and (b) 
S = 0.5 for the complete wetting case. The critical nucleus at 
the saddle point is a solid core surrounded by a liquid wetting 
layer (S = 0.5). The parameters used are a = 0.1, /? = 0.7, 
and r = 0.2. The red solid lines indicate the MFEP of the 
nucleation process. The saddle points are indicated by the 
red point on each path. 

Figures Eta) and (b) show contour plots of the free- 
energy landscape of g{x,t) in the x — t plane when the 
relative supersaturation 5 of the liquid phase is relatively 
low ((a) 6 = 0.3 and (b) S = 0.5) when a = 0.1 and 
j3 = 0.7, which satisfy the complete wetting condition 
(S/o-sv = l-(a + /8) = l- (0.1 + 0.7) = 0.2 > 0). The 
red solid line and the red point on each line indicate the 
MFEP and the saddle point, respectively. When 8 = 0.3, 
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the saddle point is still on the solid axis. However, the 
saddle point shifts from the solid axis to the liquid side 
at (x* ,t*) = (1.00, 0.92) with a free-energy barrier of nu- 
cleation g* = 0.98, which is lower than g^ v = 1.00 when 
S = 0.5. These MFEPs (red solid curves) were obtained 
by solving the overdamped equation of motion of evolu- 
tion similar to the phase-field equation^. More sophisti- 
cated numerical methods such as the string method^ are 
unsutable because the basin of attractor that corresponds 
to the bulk solid phase is located at infinity (x = oo). 

In this case (6 = 0.5), the critical nucleus at the saddle 
point is mostly a solid core surrounded by a thin layer 
of liquid. The radius of the solid core occupies 92% of 
the total radius of the composite critical nucleus at the 
saddle point. After crossing the barrier, the MFEP ap- 
proaches the bulk solid phase at (x — > oo, t — 1), and the 
composite supercritical nucleus becomes a solid nucleus. 
Therefore, the nuclcation process is not a two-step pro- 
cessor— but a one-step process with a single activation 
energy even though a macroscopically thick layer appears 
around the solid core during the evolution. 



energy barriers between the two valleys along the liquid 
axis and near the solid axis as shown in Fig. Eta). 




0.5 



1.5 



FIG. 6: Contour plot of the free-energy landscape of g(x,t) 
when the relative supersaturation is moderate (S = 0.63) for 
the complete wetting case. The parameters used are a — 0.1, 
/3 — 0.7, and r = 0.2. The two red solid curves indicate the 
two transition paths of the nucleation process. The saddle 
points are indicated by red points on the paths. There are 
two critical nuclei. One is a solid core surrounded by a liquid 
wetting layer and the other is a pure liquid nucleus. 

As we increase S further, two saddle points again ap- 
pear at (x,t) = (1.003,0.703) and (x,t) = (1.109,0.000) 
when 6 = 0.63 (Fig.[S|). The free-energy barrier at these 
two points is exactly the same (/* = 0.867). Therefore, 
two nucleation routes coexist. This situation is similar 
to that in Fig.[3][a). However, since the liquid phase can 
wet the solid nucleus, one of the critical nuclei is a com- 
posite nucleus with a solid core surrounded by a liquid 
layer whose thickness decreases as the nucleus grows after 
crossing the barrier. The other is a pure liquid nucleus. 
These two types of nucleus cannot transform between 
each other freely around the saddle point as there are 
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(b)Free-energy landscape of supercritical 
nucleus (5 = 0.63) 



FIG. 7: (a) Three-dimensional view of the free-energy surface 
shown in Fig. ([6]) near the saddle point when the relative 
supersaturation is moderate (5=0.63). Clearly there is an 
energy barrier between the valley along the liquid axis (t = 0) 
and the valley near the solid axis (t ~ 1). (b) Contour plot 
of the free-energy landscape of g(x,t) for the supercritical 
nucleus. Two MFEPs eventually converge to the upper right 
corner of the pure solid phase. 

As there is an energy barrier between the two valleys 
along t = and t ~ 1 near the saddle points, there is a 
ridge between the two valleys in Figs. [6] and [7] Thus, the 
nucleation of the stable solid from the pure liquid nucleus 
at the saddle point along the liquid axis occurs only by 
overcoming the free-energy barrier from t = to t ~ 1 
after crossing the saddle point at (x,t) = (1.109,0.000). 
Thus, the nucleation of the stable solid via the liquid 
critical nucleus appears to occur via two-step nucleation. 
In fact, the second barrier is not exactly the saddle 
point as assumed in simplified theories of two-step nu- 
clcatior>i£~— , but is the ridge. In this case, the liquid 
supercritical nucleus can survive and become long-lived 
and macroscopically large even though it is thermody- 



8 



namically mctastablc. Of course, the MFEP along this 
liquid axis should eventually merge into the MFEP along 
the solid axis far from the two saddle points because the 
free energy of the solid phase at (x — > oo, t = 1) should 
be lower than that of the liquid phase at (x — > oo, t = 0). 

Figure [T^b) shows a contour plot of the free-energy 
landscape of the supercritical nucleus far from the saddle 
point. One can easily imagine that the MFEP along the 
liquid axis eventually merges with the MFEP along the 
solid axis similarly to in Fig. SJb). In this case, however, 
we did not calculate the MFEP owing to the numeri- 
cal difficulty as the length scale is longer than that in 
Fig. HJb). Again, the nucleation process is not a two- 
step process^— but two parallel channels exist, both of 
which correspond to one-step nucleation with single ac- 
tivation energies. Thus, the nucleation rate is given by 
the formula for the parallel conductance in Eq. (|27|). 

As S is further increased and the metastable liquid 
phase becomes more stable, the saddle point shifts along 
the liquid axis to (x,t) = (0.997,0.00). Figure H] shows a 
contour plot of the free-energy landscape when S = 0.70. 
In this case, the critical nucleus is all liquid. However, 
it gradually becomes the solid after crossing the saddle 
point without overcoming the free-energy barrier. 
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FIG. 8: Contour plot of the free-energy landscape of g{x,t) 
when the relative supersaturation is high (S = 0.7) for the 
complete wetting case. The critical nucleus is all liquid along 
the liquid axis. 

As the nucleation is expected to proceed via the 
MFEP—, we can summarize the scenario of nucleation 
as follows: When the relative supersaturation of the liq- 
uid phase is low, the liquid phase is metastable and can- 
not appear, and the solid nucleus or the solid nucleus 
surrounded by a thin wetting layer of liquid directly ap- 
pears from the supersaturated vapor (Figs. (Ha), (b), and 
03a), (b)). In this case, we may expect one-step nucle- 
ation of the stable solid phase. When the supersatura- 
tion is moderate, two nucleation channels, along the liq- 
uid axis and near the solid axis, coexist (Figs. Ela) and 
I!])- When the complete wetting condition is satisfied, the 
critical nucleus near the solid axis is a composite nucleus 



that consists of a solid core surrounded by a liquid layer 
(Figs. [5] and EJa)). In this case, we may expect one- 
step nucleation again. However, we may expect that two 
nucleation processes occur in parallel and that Eq. ([27)1 
applies. Therefore, the nucleation rate is enhanced, in 
general, compared with that of two-step series nucleation 
given by Eq. (|%fy IL2g. 

When the relative supersaturation of the liquid phase 
is high, the liquid phase is relatively stable. Thus, the 
critical liquid nucleus may first appear within the super- 
saturated vapor and then pass through the saddle point 
and continue to grow (Figs. E2b), EJa), (b), and[5]). This 
growing liquid supercritical nucleus gradually changes to 
the stable solid phase far from the saddle point without 
crossing the energy barrier. This scenario of one-step 
nucleation is completely different, for example, from the 
nucleation of protein crystals ^ 15 ' 16 , where two-step nu- 
cleation with two energy barriers is assumed. 

It is well known that the observed nucleation is usually 
taking place via the metastable intermediate stated—. In 
fact, our capillarity theory is consistent to this picture 
because it predicts that the nucleation may occur via the 
metastable liquid phase when the relative supersatura- 
tion of the liquid phase is high (Figs.[3]Jb) and [8]). Since 
the metastable liquid phase is free-energetically closer to 
the stable solid phase than the metastable vapor phase, 
the nucleation occur from the vapor phase to solid phase 
via the metastable liquid phase as the liquid is almost 
stable. In this case, we may expect one step nucleation 
from vapor to liquid with a nucleation rate Jlv given by 
a single free energy barrier AG LV . Even though there 
are no true saddle point for the liquid to solid transition 
(Fig. Bib), [7|b) and |5J| after crossing the saddle point of 
the liquid nucleus and the free energy barrier in Eq. ([25)1 
appears to be zero (AGsl = 0), there will be another 
free energy barrier that comes from the preexponcntial 
factor Asl in Eq. (f25j) . Therefore the vapor to solid nu- 
cleation rate Jsv will be approximately given by Eq. ([25)) 
and will be characterized by two nucleation rates Jlv and 
Jsl- Our simple capillarity theory cannot include such 
a kinetic effect as ours is based on the quasi-equilibrium 
thermodynamics and cannot include kinetic effect. 



IV. CONCLUSION 

We have used the classical capillarity approximation 
to study the whole process of nucleation when an inter- 
mediate metastable phase is involved. By following the 
minimum-free-energy path in the free-energy landscape 
of nucleation, we studied not only the critical nucleus 
at the saddle point of the energy landscape but also the 
whole process of nucleation starting from the initial em- 
bryo. We found that the critical nucleus can be solid, 
composite, or liquid depending both on the relative su- 
persaturation of the metastable liquid phase and on the 
wetting properties of the liquid and solid phases. The 
compositions of solid and liquid in the composite nucleus 
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depends strongly on the relative supersaturation of the 
intermediate liquid phase. Although the free-energy bar- 
rier of the critical nucleus at the saddle point can also be 
studied using density functional theory, a comparison of 
the free energy only at the critical point does not reveal 
much about the nucleation process when an intermediate 
metastable phase is involved. 

Since we studied nucleation and not growth, the 
appearance of a macroscopic metastable phase during 
growth is outside the scope of this work. For such a prob- 
lem, various variants of the phase-field modelii~— will be 
useful. Finally, we stress that our simple model based on 
the capillarity approximation predicted that the inter- 
mediate metastable state (liquid) can survive even when 
the metastable liquid phase is not critica l 44 i 45 . A diver- 



gent correlation length when the intermediate metastable 
liquid phase is critica l 44 ! 45 and the long-range intermolec- 
ular interaction 4 ^ will certainly affect the conclusions de- 
rived from our simplified capillarity theory of nucleation 
based on the short-range interaction. These issues are 
left for future investigations. 
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